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ABSTRACT 



Aims. We study the visibility of sunspots and its influence on observed values of sunspot region parameters. 

Methods. We use Virtual Observatory tools provided by AstroGrid to analyse a sample of 6862 sunspot regions. By studying the 
distributions of locations where sunspots were first and last observed on the solar disk, we derive the visibility function of sunspots, 
the rate of magnetic flux emergence and the ratio between the durations of growth and decay phases of solar active regions. 
Results. We demonstrate that the visibility of small sunspots has a strong centre-to-limb variation, far larger than would be expected 
from geometrical (projection) effects. This results in a large number of young spots being invisible: 44% of new regions emerging in 
the west of the Sun go undetected. For sunspot regions that are detected, large differences exist between actual locations and times of 
flux emergence, and the apparent ones derived from sunspot data. The duration of the growth phase of solar regions has been, up to 
now, underestimated. 
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1. Introduction 

The birth of a new spot on the solar disk indicates the emergence 
of magnetic flux through the photosphere, a process which is 
key to the solar cycle (Solanki 2003; Fisher et al. 2000) and the 
study of stellar magnetic dynamos. Sunspots also cause varia- 
tions in the total solar iiTadiance, an important parameter in de- 
termining the Sun's influence on climate (Foukal et al. 2006). 
The presence of a sunspot is key to a solar region being as- 
signed an active region number by the NOAA Space Weather 
Prediction Center (http://www.swpc.noaa.gov) so that its evolu- 
tion and activity can be tracked (Gallagher et al. 2007). The for- 
mation and evolution of active regions are fundamental to solar 
dynamic phenomena such as flares and Coronal Mass Ejections 
and their effect on the Earth environment. 

The visibility of sunspots is currently thought to be limited 
only by geometrical effects arising from projection of the solar 
sphere onto a 2D image, effects referred to as foreshortening. 
In this paper we present results obtained serendipitously while 
analysing sunspot data by means of Virtual Observatory tools, 
showing that the visibility of small sunspots is much poorer than 
predicted by the foreshortening model. 



2. Data analysis 

We analysed sunspot group data from the USAF/Mount 
Wilson catalogue, for times between 1 December 1981 to 31 
December 2005, covering two and a half solar cycles. We pro- 
cessed the data by means of AstroGrid workflows. AstroGrid 
(http://www.astrogrid.org) (Walton et al. 2006) is the UK's con- 
tribution to a global Virtual Observatory ( VO), aimed at allowing 



seamless access to a variety of astronomical data, and at pro- 
viding efficient software tools for data analysis. Our workflows 
processed all entries in the catalogue and extracted properties of 
individual regions, identified by their NOAA region number, in- 
cluding the time and location of their first and last observation. 
The longitude of each region at 12:00UT on the days when it 
was first and last observed was calculated. 

Figure 1 shows histograms of the locations at which sunspot 
regions were first and last observed, with 0° longitude corre- 
sponding to the Earth-Sun line. As the Sun rotates, many spots 
first come into view near the east (or 'rising') limb, causing the 
peak to the left in Fig. 1-a. Similarly, the peak in Fig. 1-b cor- 
responds to regions that rotated out of view. Regions first ob- 
served sufficiently far from the east limb, are generally assumed 
to be 'new', indicating the emergence of magnetic flux through 
the photosphere. As spots move in longitude by approximately 
14.3° each day, regions first observed to the west of -60° are 
typically described as new emergences (see inset in Fig. 1-a). 
Similarly regions last seen to the east of -1-60° are interpreted as 
having decayed to the point that a spot is no longer visible (see 
inset in Fig. 1-b). 

The insets in Fig. 1 display strong east-west asymmetry. 
A total of 825 new regions are seen to emerge in the bin 
[-60°,-40°], while only 177 in [h-40°,-h60°], a ratio of 4.7:1. 
What is the cause of the strong asymmetry in these curves? Why 
should the number of new regions emerging in any given longi- 
tude bin not be constant? What is the true rate of magnetic flux 
emergence on the Sun? 

An asymmetry in the location of emergence of new 
sunspots as viewed from Earth was discovered 100 years ago 
(Maunder 1907) and attracted the attention of famous physicists. 
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Fig. 1. Histograms of number of sunspot regions versus their longitude 
at 12:00 UTC on the day when they were first (a) and last (b) observed . 
Each longitude bin is 6°. Data are from the USAF/Mt Wilson catalogue 
of sunspot groups, for the time range 1 December 1981 to 31 December 
2005. The total number of regions is 6862. Longitude 0° corresponds 
to the Earth-Sun line, negative values to eastern longitudes and positive 
values to western longitudes. Insets show the histograms for the range 
[-60°, +60°]. 



who demonstrated that a visibility function that favours obser- 
vations in the centre of the disk, and the curve of evolution 
of a spot's size, can produce such asymmetry (Schuster 1911; 
Minnaert 1939). The graphical representation introduced by 
Minnaert (1939) makes the cause of the asymmetry immedi- 
ately clear. However, Minnaert himself went on to assume that 
the only factor limiting the visibility of sunspots is geometrical 
effects associated with foreshortening. This results in a visibil- 
ity function proportional to sec A, with A the longitude, and poor 
visibility only very near the limb. The latter assumption has re- 
mained undisputed until the present day. Moreover, appreciation 
of the cause of the east-west asymmetry appears to have been 
lost, and asymmetries in sunspot parameters have since been 
ascribed, e.g., to observer bias or systematic inclination of the 
magnetic flux tubes (Howard 1991). 

The large asymmetries in the insets of Fig. 1, however, 
demonstrate that the visibility of new sunspot formation in 
[-60°, -1-60°] is poor, and flux appears to emerge and disperse at 
locations (and times) that are far from the locations (and times) 
where actual emergence and decay took place. 

We use the data of Fig. 1 to derive the rate of region emer- 
gence, the average ratio between the slopes of growth and decay 
phases of region evolution and to determine the visibility func- 
tion. Let s(A) be the visibility function, giving the minimum ac- 
tual (as opposed to apparent) area that a sunspot region needs to 
reach to be visible at longitude A. s is expected to have a mini- 
mum at A-Q, where visibility is best. Assuming that, over a given 
time period, the number of magnetic flux emergences in a unit 
longitude bin is a constant A^i , the number N{A) of regions ob- 
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Fig. 2. Values of the sums of counts in positive and negative longitude 
bins at either side of the Earth-Sun line, as a function of absolute value 
of the longitude, for region emergences (a) and disappearances (b). 
Dotted lines indicate the mean. The sizes of error bars are also shown. 
The size of each bin is 6°. 



served emerging in a unit bin at longitude A is (Schuster 1911) 
[See also Appendix A]: 



N(A) = A^i 



1 



s'(A) 



(1) 



where Q. is the solar rotation rate and k is the linear growth con- 
stant of a region's area with time. Similarly, we obtain for the 
number n(A) of regions seen to disappear at A: 



n(A) — «! 



Q 

1 + - s'iA) 



(2) 



with «! the number of regions that reach peak area in a unit lon- 
gitude interval and / the decay time constant, Z>0. 

Considering two bins centred at +A and -A and indicating 
the number of regions seen emerging in them as and re- 
spectively, one can write Eq. (1) for each of the two bins, and 
obtain expressions for and A^_. By adding these together and 
assuming that the visibility function is symmetric with respect 
to A-0, so that s'(-A)--s'(+A), one finds that N++N--2Ni 
(Schuster 1911) [a misprint appears in Schuster's expression on 
his p. 319]. For regions disappearing in the same two bins, from 
Eq.(2): n+ + n--2ni-2Ni. Here we assume that the average 
number n\ of regions that reach peak area in a unit longitude 
bin is equal to the average number A^i of regions that emerge in 
a unit longitude bin. 

Figure 2-a shows values of N++N- versus \A\ obtained from 
the data of Fig. 1-a by adding counts in positive and negative 
/l-bins at either side of A-Q. Figure 2-b shows n++n^ obtained 
from the disappearences data of Fig. 1-b. N++N- and n+-i-n_ are 
approximately constant, as predicted by Schuster's theory, and 
we can use their mean values (indicated by the dotted lines in 
Fig. 2) to calculate A^i, the actual number of sunspot regions 
emerging in a unit bin. We find A^i = 160.55+1 1.41 from Fig. 2-a 
and 158.95+12.19 from Fig. 2-b. Considering that all regions 
in the USAF/Mt Wilson catalogue are included in our analy- 
sis, and that the catalogue is compiled from observations made 
at discrete intervals of time, the data fit Schuster's simple the- 
ory remarkably well. There is excellent agreement in the values 
of A^i derived separately from emergences and disappearances. 
The data shown in Fig. 2 are consistent with a constant value 
of N++N- and n++n- as would be expected from a visibility 
function that is symmetric with respect to A=0. Figure 2-a does 
display a small slope, while this is not seen in Fig. 2-b. Several 
authors discussed the issue of whether the magnetic flux tubes 
of active regions present a systematic inclination with respect to 
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the direction perpendicular to the solar surface. An analysis of 
magnetograms showed a systematic inclination of magnetic flux 
tubes of growing active regions of 24° in the W-E direction, i.e. 
trailing the rotation (Howard 1991). This would result in better 
visibility of young regions in the west, the opposite of the effect 
we find. We conclude that our data do not show evidence for 
a strong inclination although this issue may need to be further 
investigated with other data. 

From Eqs.(l) and (2), by solving for the first derivative s'(A), 
and by dividing one equation by the other, we find an expres- 
sion for the ratio k/l between the growth and decay constants 
characterising sunspot region evolution. We find a mean value 
fc//=l. 37+0.26 (s.d.), from 16 longitude bins. (Here, we exclude 
4 longitude bins near A=0 in which N/N] and n/n\ are close to 
1, giving two terms close to zero to be divided by each other to 
find k/l). Our value for k/l implies that, on average, the dura- 
tion of the decay phase of an active region is only 1.37 times 
larger than that of the rise phase. This is very different from pre- 
vious estimates, according to which active regions reach maxi- 
mum development very quickly, e.g. within 5 days of emergence 
(Harvey 1993), and decay slowly. The ratio being different from 
1 is the cause of the lack of complete symmetry in the distribu- 
tions of emergences and disappearances (see insets of Fig. 1). 
While in general the curve of evolution of sunspot regions wiU 
be complex, its hnear envelope, describing a triangle with slopes 
k and /, is a useful first approximation. By means of modelUng, 
we found that the asymmetry in [-60°, +60°] is largely indepen- 
dent of the actual total lifetime of regions, while lifetime is the 
key factor influencing the height of the peak in Fig. 1-a. 

Having obtained A^i and the ratio k/l, values of the deriva- 
tive of the visibility function, s', given by Eqs.(l) and (2), de- 
pend on a single parameter, the growth constant k. In Fig. 3-a 
we plot s' data points, calculated using Eqs.(l) (squares) and 
(2) (triangles) for a value of the growth constant k-40 msh/day 
and using k/ 1=131 (msh indicates milhonths of the visible so- 
lar hemisphere). We fit the function s'(A)= c\ tan"' (c2 X) to the 
data and obtain a best fit when ci =117.0 and C2=4.7 (the solid 
line in Fig. 3-a). By integration we obtain s(A)=ci /c2 [x tan ' x- 
1 /2 ln(l + x^)] + A„,„, where x=C2 A (with A in radians) and A„,„ 
is the minimum of the visibility function, set equal to 8 msh, 
the approximate minimum area required at the disk centre for a 
group to be included in the catalogue. Figure 3-b shows s(A) for 
the fit of Fig. 3-a (solid fine) and for those obtained for two other 
values of k. The dotted line shows A^m sec A, the visibiUty func- 
tion expected from geometrical considerations describing fore- 
shortening. 

Sunspot regions are in reality characterised by a distribution 
of growth and decay constants. The curves of Fig. 3-b therefore 
need to be interpreted as the visibility functions arising when 
the most probable value of k in the distribution is equal to the 
numerical value given. Recurrent sunspots were found to have 
a decay rate constant 1^10 msh/day (Martinez Pillet et al. 1993), 
giving k=l4 msh/day using our k/l ratio. Due to the very slow 
decay of recurrent spots the latter is a hkely good lower limit for / 
and consequently k. Even when fc= 14 msh/day, sunspot visibility 
is much worse than expected from projection effects only. To 
recover the geometrical visibility curve from the data, it would 
be necessary to assume k=3 msh/day, a value not consistent with 
observations. 

3. Discussion 

Figure 3-b demonstrates that the visibility of small sunspots is 
much poorer than expected from geometrical effects associated 
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Fig. 3. (a) Derivative of the visibility function from emergence (squares) 
and disappearance (triangles) data, for k=40 msh/day and k/l=\31, with 
msh=millionths of the visible solar hemisphere; the solid line shows 
the best fit to the data points, (b) Visibility function s(X), giving the 
minimum area a sunspot needs to reach to be detected at longitude X, 
for several values of k. The dotted line shows the function A„i„ sec X. 



with foreshortening. It shows that the minimum area required for 
a spot to be detected at /l=±30° is more than twice the thresh- 
old area at A={f if A:=14 msh/day and almost 4 times if A:=30 
msh/day. The centre-to-limb variation in visibility is remarkably 
large. 

We investigated whether the asynnmetries shown in Fig. 1 
display any solar cycle dependence, and found no evidence of 
it. The USAF/Mount Wilson locations and times of first appear- 
ances agree with those from SOHO/MDI continuum data (as ver- 
ified manually for a sample of regions). It is known that seeing 
associated with ground-based data does not cause significant re- 
duction of visibility (Gyori et al. 2004). 

How many spots are affected by the visibility effects here 
described? The distribution of sunspot areas measured at a 
single longitudinal location on the solar disk is lognormal 
(Bogdan et al. 1988) and the number of spots with area at 
Central Meridian around 10 msh is more than 2 orders of magni- 
tude larger than the number of spots with area of about 100 msh. 
Therefore the majority of sunspots wiU cross the visibility curve 
shown in Fig. 3-b. 

Our results have a number of important implications. The 
first is that the radiative processes that make a small region 
of strong magnetic field appear as a dark spot, have a strong 



4 



Dalla et al.: Invisible sunspots 



centre-to-limb variation. This may prove important for the study 
of sunspots' 3D structure and will require further investiga- 
tion. Whether larger spots are also affected by the same pro- 
cess will also need further study. Reports of centre-to-limb 
variations of corrected sunspot areas have appeared in the 
literature (Gyori et al. 2004; Hoyt et al. 1983). Faculae have a 
large centre-to-limb variation, and their contrast changes sign 
as one moves towards the disk centre, resulting in their be- 
ing darker than the surrounding photosphere at the disk centre 
(Lawrence et al. 1993). This demonstrates that the appearance of 
photospheric magnetic flux tubes strongly depends on the view- 
ing angle. 

The second implication is that actual distributions of 
sunspot lifetimes and areas may differ from the ap- 
parent ones derived from observations. The latter have 
been used to constrain mechanisms of sunspot formation 
and decay (Solanki 2003; Petrovay & Van Driel-Gesztelyi 1997; 
Martinez Pillet et al. 1993). A large number of regions reported 
of short duration may in fact have longer lifetimes, and be cross- 
ing in and out of the visibility curve. The Gnevyshev-Waldmeier 
law, stating that a sunspot's lifetime increases linearly with its 
maximum size, may need to be reassessed in light of our results. 
Poor visibility of region emergence means that the actual time of 
magnetic flux emergence can be much earlier than the apparent 
time, e.g. for a region seen to emerge at /l=-50°, by approxi- 
mately 2 days. On the other hand, a large fraction of new emer- 
gences in the western portion of the solar disk go undetected. By 
using the data in the inset of Fig. 1-a for A>Q and the value of 
actual number of emergences A^i=160 obtained from the data of 
Fig. 2, we obtain that 44% of new spots emerging in [0°, -h60°] 
were invisible. 

The presence of a sunspot is key to a solar active region be- 
ing assigned an Active Region Number by the NOAA Space 
Weather Prediction Center (see e.g. Dalla et al. (2007) for the 
full list of criteria). A region that has been given a NOAA num- 
ber is monitored and its activity tracked (Gallagher et al. 2007). 
New regions emerging in the west of the Sun with their spots be- 
ing invisible are missed and not tracked. We conclude that cur- 
rent criteria for assigning Active Region Numbers may need to 
be revised and that EUV solar images may need to be routinely 
used to supplement white-light information. 

Sunspots are well known to cause depletions in the total 
solar irradiance (TSI) (Foukal et al. 2006) and many models of 
TSI variation need as input information on the number and ar- 
eas of sunspots. While the sunspots that most affect TSI are the 
largest ones, of area typically well above the visibility thresh- 
old shown in Fig. 3-b, our results impact TSI studies because 
they demonstrate that the apparent age and stage of development 
of a sunspot may not correspond to the actual ones. The latter 
information is required when studying the time dependence of 
sunspot effects on TSI and whether the age of a region is an im- 
portant factor in determining the magnitude of TSI decrease. 

The asymmetry in the distribution of emergences was ob- 
tained, unexpectedly, during a study aiming at cross correlating 
catalogues of sunspot regions and flares, by means of AstroGrid 
workflows (Dalla et al. 2007). This demonstrates the usefulness 
of VO tools in making new science possible, by provision of 
better tools for analysis of large datasets. 
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Appendix A: Schuster's equation and l\/linnaert's 
graphical representation 

The fact that a visibility function favouring sunspot observa- 
tions in the centre of the solar disk, should result in an east- 
west asymmetry in the number of regions seen emerging at the 
Sun, is not immediately intuitive. In this Appendix we sum- 
marise Schuster's derivation of Eq.(l) of our paper (Schuster 
1911) and describe Minnaert's graphical representation, from 
which the cause of the asymmetry becomes immediately appar- 
ent (Minnaert 1939). 

Two phenomena combine to produce the effect here de- 
scribed. The first is the fact that sunspot regions are evolving: 
their evolution can be characterised in terms of their area and in 
a zero-th order approximation can be described by a curve such 
as shown in Fig. A.l: a growth phase with slope k and a decay 
with slope -/ (Z>0). The second phenomenon is the solar rota- 
tion, which carries sunspots which emerged in the east of the 
Sun towards regions of better visibility (i.e. towards the centre 
of the disk) and regions that emerge in the west towards regions 
of worse visibility. 

The combination of region evolution and rotation, together 
with the specific form of the visibility function, determines 
whether a sunspot region is seen or not, and the location and 
time of its first appearance (and of its disappearance) to an Earth 
observer. This is clear from Minnaert's graphical representation 
(Minnaert 1939), as shown in Fig. A. 2. The curve in the dia- 
gram represents the visibility function s{A), giving the actual 
area that a sunspot region needs to reach to be visible at lon- 
gitude A. The parallel lines at a slope represent growth phases 
of sunspots. When a sunspot region growing along a given line 
crosses s{X), it becomes visible. From this graphical represen- 
tation it can be seen that if the visibility function has a strong 
centre-to-limb variation, many spots forming in the west are in- 
visible. An east-west asymmetry in the number of regions ob- 
served emerging thus results, and the asymmetry depends on the 
gradient s' of the visibility function. 

The process qualitatively represented by Minnaert's graph 
was quantitatively described by Schuster (1911), as follows. Let 
Al indicate the longitude at which a sunspot forms, and A the 
longitude at which it is first seen because it has crossed the visi- 
bility curve s{A). The two longitudes are related by the following 
equation: 



A-A^ 



^s{A) 



(A.l) 



where Aq is the area of the spot when it first emerges, k is the 
slope of the growth phase and Q the solar rotation rate. Eq. A. 1 




Longitude 

Fig. A.2. Schematic of Minnaert's graphical representation (adapted 
from Minnnaert 1939). The black curve is the visibility function s{X) 
and the grey parallel curves represent the growth curves of sunspot re- 
gions. 



expresses the fact that the area at the location where the spot is 
first seen is equal to s{A). 

Eq. A. 1 can be re-arranged to give: 



Q Q 

A- — s{A) ^Ai-—Aq 

which, differentiated, gives: 



dA 



Q , 
1 - - s'(A) 
k 



dA] 



(A.2) 



(A.3) 



If A^i indicates the number of sunspots emerging in a unit 
longitude bin and N{ A) the number of regions observed emerging 
in a unit bin at longitude A, then A^i dA\ -N{A) dA, which gives, 
using Eq. A.3: 



Q. , 
1 - -- s\A) 
k 



(A.4) 



(Schuster 1911) (Eq (1) of our paper). 

Applying the same derivation to the decay phase, charac- 
terised by a slope -/, we obtained Eq. (2) of our paper 




Fig. A.l. Schematic of time evolution of a sunspot's area. Here Aq is 
the area of the spot at the time when it first emerges through the photo- 
sphere. 



